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ABSTRACT 



We have succeeded in obtaining a precision analytic fitting formula for the exact numerical 
results of the relativistic corrections to the thermal Sunyaev-Zeldovich effect for clusters of 
galaxies which has a 1% accuracy for the crossover frequency region where the thermal thermal 
Sunyaev-Zeldovich effect signal changes from negative to positive sign. The fitting has been 
carried out for the ranges 0.020 < e < 0.035 and < X < 15, where 8 e = kBT e /m e c 2 , 
X = huj/ksTo, T e is the electron temperature, u> is the angular frequency of the photon, and 
To is the temperature of the cosmic microwave background radiation. The overall accuracy of 
the fitting is better than 0.1%. The present analytic fitting formula will be useful for accurate 
analyses of the thermal Sunyaev-Zeldovich effect for clusters of galaxies. 



Subject headings: cosmic microwave background — cosmology: theory — galaxies: clusters: 
general — radiation mechanisms: thermal — relativity 
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1. INTRODUCTION 

In a recent series of papers the present authors (Itoh, Kohyama, & Nozawa 1998; Nozawa, Itoh, & 
Kohyama 1998b; Itoh, Nozawa, & Kohyama 2000a) have calculated the relativistic corrrections to the 
thermal, kinematical, and polarization Sunyaev-Zeldovich effects (Zeldovich & Sunyaev 1969; Sunyaev & 
Zeldovich 1972, 1980a, 1980b, 1981). Other recent references on this subject include Rcphaeli (1995); 
Rephaeli & Yankovitch (1997), Stebbins (1997), Challinor & Lasenby (1998,1999), Sazonov & Sunyaev 
(1998a, 1998b, 1999). The present authors have also calculated the relativistic thermal bremsstrahlung 
Gaunt factor for the intercluster plasma (Nozawa, Itoh, & Kohyama 1998a; Itoh et al. 2000b). 

The present authors (Nozawa et al. 2000) have presented an accurate analytic fitting formula for the 
fractional distortion of the photon distribution function caused by the thermal Sunyaev-Zeldovich effect 
which has been calculated by the direct numerical integration of the collision term of the Boltzmann 
equation (Itoh, Kohyama, & Nozawa 1998). They have obtained 0.1% accuracy of the analytic fitting 
formula for most of the parameter region 0.02 < e < 0.05 and < X < 20, where 9 e = kBT e /m e c 2 , 
X = Tiw/fcsTo, T e is the electron temperature, w is the angular frequency of the photon, and T is the 
temperature of the cosmic microwave background radiation. However, the relative error in the crossover 
frequency region where X « 3.9 amounted in some cases to more than 5% owing to the fact that the 
fractional distortion of the photon distribution function caused by the thermal Sunyaev-Zeldovich effect is 
such a rapidly varying function of X in this region that the excellent overall fitting is extremely difficult. 

In view of the recent spectacular advances in the accuracy of the observation of the Sunyaev-Zeldovich 
effect for clusters of galaxies (Birkinshaw 1999), it is worthwhile to present a still more accurate fitting 
formula for the thermal Sunyaev-Zeldovich effect which has an overall 0.1% accuracy and has a 1% accuracy 
for the crossover frequency region X « 3.9. In order to achieve this goal we will restrict the range of the 
electron temperature to 9 e = kBT e /m e c 2 < 0.035. This temperature range covers all the observed galaxy 
clusters. 

The present paper is organized as follows. In § 2 we give the the analytic fitting formula. Concluding 
remarks will be given in § 3. 
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2. ANALYTIC FITTING FORMULA 

Itoh, Kohyama, & Nozawa (1998) have calculated the fractional distortion of the photon distribution 
function caused by the thermal Sunyaev-Zeldovich effect An(X)/n n (X) by numerical integration of the 
collision term of the Boltzmann equation: 



An(X) 
n (X) 



yF(9 e ,X), (2.1) 



y = <t t J d£N e , (2.3) 

Oe = ^f, (2.4) 

where u> is the angular frequency of the photon, and T is the temperature of the cosmic microwave 
background radiation, ot is the Thomson scattering cross section, N e is the electron number density, T e is 
the electron temperature, and the integral in equation (2.3) is over the path length of the galaxy cluster. 
Nozawa et al. (2000) have presented an analytic fitting formula for F(6 e ,X). Since F(6 e ,X) is a rapidly 
varying function of X, it was extremely difficult to obtain an excellent overall fitting. Nozawa et al.'s (2000) 
fitting formula has a general accuracy of 0.1%. However, in the crossover frequency region where X w 3.9, 
the relative error in some cases amounted to more than 5%. In view of the recent spectacular advances in 
the accuracy of the observation of the Sunyaev-Zeldovich effect for clusters of galaxies (Birkinshaw 1999), 
it is worthwhile to present a still more accurate fitting formula. 

In the observation of the thermal Sunyaev-Zeldovich effect, the directly observed quantity corresponds 
to the distortion of the spectral intensity 

X 3 An{X) 



AI = 



e x - 1 n (X) 



= y^jF(e e ,X). (2.5) 
In this paper we will present an analytic fitting formula for A J instead of An(X)/no(X). 

We express the fitting formula for 0.020 < 9 e < 0.035, < X < 15 as follows: 

AI 1 X 3 An(X) 

y y e x - 1 n (X) 

X 3 Xe x 

= -ir-^ ^x~^ ( y o + o e Yi + e 2 e Y 2 + e 3 e Y 3 + e 4 e Y i ) + r. (2.6) 

e — 1 e — 1 

The functions Y , Y\, Y 2 , Y 3 , and Y4 have been obtained by Itoh, Kohyama, and Nozawa (1998) with the 
Fokkcr-Planck expansion method, and their explicit expressions have been given. The graph of the redidual 
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function R is shown in Figures 1,2. We take the fitting formula for R as 

, for < X < 2.4 



R = { ™ (2.7) 
' a l3 &l Z 1 , for 2.4 < X < 15.0 , 

i,j=0 



where 



e e = ^(6> e -o.o2), 0.020 < e e < 0.035, (2.8) 

Z = JL(X-8.7), 2.4<X<15.0. (2.9) 

The coefficients dij are presented in Table 1. The accuracy of the fitting formula (2.6), (2.7) is generally 
better than 0.1%. The fitting formula has 1% accuracy in the crossover frequency region where X « 3.9. 
For 6 e < 0.02, the results of Itoh, Kohyama, & Nozawa (1998) give sufficiently accurate results (the accuracy 
is generally better than 1% for < X < 15). 



3. CONCLUDING REMARKS 

We have presented a precision analytic fitting formula for the spectral intensity distortion caused by 
the thermal Sunyaev-Zeldovich effect which has been calculated by the direct numerical integration of 
the collision term of the Boltzmann equation (Itoh, Kohyama, & Nozawa 1998). The spectral intensity 
distortion corresponds to the directly observed quantity. The fitting has been carried out for the ranges 
0.020 <0 e < 0.035, < X < 15. The accuracy of the fitting is generally better than 0.1% and is about 1% 
in the crossover frequency region where X « 3.9. The present results will be useful for accurate analyses 
of the thermal Sunyaev-Zeldovich effect for clusters of galaxies. For galaxy clusters with relatively low 
temperatures 9 e < 0.02, the Fokker-Plank expansion results of Itoh, Kohyama, & Nozawa (1998) will be 
sufficiently accurate (the accuracy is generally better than 1%). 

We thank Professor Y. Oyanagi for allowing us to use the least-squares fitting program SALS. This 
work is financially supported in part by a Grant-in- Aid of Japanese Ministry of Education, Science, Sports, 
and Culture under contract 10640289. 
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Figure Captions 

• Fig.l. The residual function R in equation (2.6) as a function of X for 9 e = 0.02. 

• Fig.2. Same as Fig.l, but for 9 e = 0.03. 
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Table 1 Coefficients a 
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i = 3 




i = 4 


i = 5 









2.38348E- 
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-4 


-4.87242E- 
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Table 1 (Continued) Coefficients a 
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2 


1.18551E- 
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